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ABSTRACT 

An  adaptive  radio-frequency  hyperthermia  system  for  improved  therapeutic 
tumor  heating  is  investigated.  Adaptive  array  techniques  are  used  to  modify  the 
electric-field  and  temperature  distribution  in  hyperthermia  experiments  with  phan¬ 
tom  targets.  A  commercial  hyperthermia  phased-array  antenna  system  at  the  SUNY 
Health  Science  Center  in  Syracuse,  New  York,  has  been  modified  to  implement 
adaptive  nulling  and  adaptive  focusing  algorithms.  The  hyperthermia  system  is  the 
BSD  Medical  Corporation  Model  BSD-2000  with  Sigma-60  annular  phased-array 
antenna  applicator.  The  applicator  operates  from  60  to  120  MHz  and  consists  of 
four  pairs  of  dipole  antenna  radiators.  The  four-channel  transmit  array  is  made 
adaptive  by  software  modifications  which  invoke  a  gradient-search  feedback  algo¬ 
rithm.  The  gradient-search  algorithm  implements  the  method  of  steepest  descent 
for  adaptive  nulling  and  the  method  of  steepest  ascent  for  adaptive  focusing.  The 
feedback  signals  are  provided  by  electric-field  short-dipole  probe  antennas. 

With  the  proposed  adaptive  hyperthermia  array  design  concept,  it  may  be 
possible  to  maximize  the  applied  electric  field  at  a  tumor  position  in  the  target  body 
and  simultaneously  minimize  or  reduce  the  electric  field  at  target  positions  where 
undesired  high-temperature  regions  (hot  spots)  occur.  The  measured  phantom- 
target  data  indicate  that  adaptive  nulling  can  reduce  the  electric  field  at  one  or  more 
target  positions  while  simultaneously  focusing  at  a  deep-seated  position  within  the 
target. 
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1.  INTRODUCTION 


Adaptive  array  antennas  are  well  known  for  their  ability  to  improve  the  performance  of 
communications  and  radar  systems  [1,  2].  Recently,  adaptive  array  techniques  have  been  applied  to 
the  medical  application  of  radio-frequency  (RF)  hyperthermia  for  deep-seated  tumor  therapy  [3-6]. 
With  an  adaptive  radio-frequency  hyperthermia  array,  it  is  theoretically  possible  to  automatically 
control  the  RF  electric  field  at  multiple  positions  within  the  target  body  [3-6].  The  electric  field 
radiated  hy  a  hyperthermia  phased  array  can  be  minimized  (nulled)  and  m2Lximized  (focused) 
at  desired  tuget  positions  by  adaptively  adjusting  the  transmit  amplifiers  and  phase  shifters  of 
the  hyperthermia  apparatus.  In  this  report,  multiple  adaptive  electric-field  nulls  and  adaptive 
focusing  are  experimentally  demonstrated  at  the  State  University  of  New  York  (SUNY)  Health 
Science  Center  with  a  modified  commercial  hyperthermia  phased-array  system  and  phantom  target 
measurements  [7]. 

The  treatment  of  deep-seated  malignant  tumors  within  a  patient  is  often  a  difficult  task.  The 
objective  of  the  treatment  is  to  reduce  in  size  or  completely  remove  the  tumor  mass  by  one  or  more 
modalities  available  at  the  treatment  facility.  Common  modalities  are  surgery,  chemotherapy,  zmd  x- 
ray  therapy  [8].  A  particular  modality  used  alone  or  in  conjunction  with  another  modality  is  “tissue 
heating,”  or  hyperthermia  [8-12].  Hyperthermia  can  be  considered  as  a  form  of  high  fever  within 
the  body;  a  controlled  thermal  dose  distribution  is  required  for  hyperthermia  to  have  a  therapeutic 
value.  Typical  localized-hyperthermia  temperatures  required  for  therapeutic  treatment  of  cancer 
are  in  the  42.5°-45®C  range.  Normal  tissue  should  be  kept  at  temperatures  below  42.5°C  during  the 
treatment.  The  most  difficult  aspect  of  implementing  hyperthermia,  with  either  radio-frequency 
waves  or  acoustic  (ultrasound)  waves,  is  producing  sufficient  heating  at  depth.  Multiple-applicator 
RF  hyperthermia  arrays  are  commonly  used  to  provide  a  focused  main  beam  at  the  tumor  position. 
A  focal  region  should  be  concentrated  at  the  tumor  with  minimal  energy  delivered  to  surrounding 
normal  tissue.  As  the  hyperthermia  antenna  beamwidth  is  proportional  to  the  wavelength,  a  small 
focal  region  suggests  that  the  RF  wavelength  be  as  small  as  possible.  However,  due  to  propagation 
losses  in  tissue,  the  RF  depth  of  penetration  decreases  with  increasing  transmit  frequency.  One  of 
the  major  problems  in  heating  a  deep-seated  tumor  with  a  hyperthermia  antenna  is  the  formation  of 
undesired  “hot  spots”  in  surrounding  tissue.  This  additional  undesired  heating  often  produces  pain, 
burns,  and  blistering  in  the  patient,  which  requires  terminating  the  treatment.  Thus,  techniques 
for  reducing  hot  spots  are  necessary  in  hyperthermia  treatment. 

A  typical  clinical  RF  hyperthermia  treatment  consists  of  several  two-hour  sessions  spread 
over  a  period  of  several  weeks.  During  a  two-hour  session,  the  first  hour  is  used  for  patient  prepa¬ 
ration,  during  which  temperature  and  vital-signs  monitoring  instrumentation  are  interfaced  with 
the  patient.  After  the  patient  is  readied,  the  hyperthermia  antenna  RF  signal  is  turned  on  and 
the  tumor  is  heated  by  the  radiated  electromagnetic  energy.  The  tumor  should  be  heated  for  ap¬ 
proximately  45  minutes  at  a  temperature  of  42.5*’-45®C.  Approximately  15  minutes  (maximum)  is 
available  prior  to  the  actual  delivery  of  RF  energy  to  the  tumor  for  manually  or  adaptively  shaping 
the  RF  electric-field  distribution  within  the  patient  target  tissue.  Current  clinical  operation  of 
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hyperthermia  phased  arrays  allows  limited  manual  control  of  the  array  transmit-element  gain  and 
phase.  Some  improvement  in  the  electric-field  distribution  can  be  achieved  by  this  manual  trial- 
and-error  method,  but  automatic  adjustment  techniques,  offered  by  adaptive  arrays,  are  desirable 
and  possibly  offer  better  electric-field  distributions. 

Several  journals  have  published  special  issues  on  acoustic  and  electromagnetic  hyperther¬ 
mia  treatment  of  cancer  [13,  14].  Electromagnetic  transmitting  arrays  in  the  frequency  band  60- 
2000  MHz  are  used  to  localize  heating  of  malignant  tumors  within  a  target  body.  Many  studies 
have  been  conducted  to  produce  improved  therapeutic  field  distributions  with  hyperthermia  phased 
arrays  [15-32].  Phase  control  can  be  used  to  synthesize  improved  RF  radiation  patterns  without 
adaptive  control  of  the  transmit-array  weights  [15-20].  Array  transmit  weights  can  be  adaptively 
controlled  to  maximize  the  tumor  temperature  (or  RF  power  delivered  to  the  tumor)  while  minimiz¬ 
ing  the  surrounding  tissue  temperature  (or  RF  power  delivered  to  the  surrounding  tissue)  [21-32]. 
All  other  studies  require  invasive  techniques  to  optimize  the  radiation  pattern  [21-32].  A  previous 
report  [3]  and  several  conference  articles  [4-6]  investigate  the  theoretical  benefit  of  using  adaptive 
nulling  with  noninvasive  auxiliary  dipole  sensors  to  reduce  the  field  intensity  at  selected  positions 
in  the  target  body  while  maintaining  a  desired  focus  at  a  tumor.  Multiple  adaptive  electric-field 
nulls  were  used  to  show  a  theoretical  reduction  in  hot  spots  for  a  homogeneous  elliptical  phantom 
target  surrounded  by  a  water  bolus  and  hyperthermia  ring  array  [3-6]. 

The  hyperthermia  annular  phased-array  antenna  system  (Model  BSD-2000,  SlGMA-60  ap¬ 
plicator  [33],  BSD  Medical  Corporation,  Salt  Lake  City,  Utah)  used  in  these  experiments  is  shown 
in  Figure  1.  *  By  fully  surrounding  the  patient  with  an  annular  phased  array,  it  is  possible  to 
obtain  constructive  interference  (or  signal  enhancement)  deep  within  the  target  volume.  The  BSD- 
2000  hyperthermia  system  uses  a  59-cm  ring-array  diameter  with  eight  uniformly  spaced  dipole 
elements  operating  over  the  frequency  band  60-120  MHz  [34].  The  eight  dipoles  are  fed  as  four 
active  pairs  of  elements.  There  are  four  high-power  amplifiers  which  drive  the  dipole  pairs  with 
up  to  500  W  average  power  per  channel.  Each  of  the  four  active  channels  has  an  electronically 
controlled  variable-phase  shifter  for  focusing  the  array.  Temperature  and  electric-field  probe  sen¬ 
sors  (both  invasive  and  noninvasive)  are  used  to  monitor  the  treatment.  A  cool-water  (5°-40®C) 
bolus  between  the  patient  and  the  phased  array  is  used  to  prevent  excess  heating  of  the  skin  surface 
during  the  treatment.  The  bolus  is  filled  with  circulating  deionized  water,  which  has  a  very  low 
RF  propagation  loss. 

A  candidate  adaptive-nulling  algorithm  is  a  gradient  search  based  on  minimizing  the  signal  re¬ 
ceived  by  electric-field  sensors  at  the  desired  null  petitions  while  maintaining  the  signal  amplitude  at 
the  foc2d  point.  The  reason  for  considering  this  algorithm  is  that  the  BSD- 2000  system  measures  the 
applied  electric-field  amplitude  but  not  the  phase.  Note;  the  E-field  probes  used  in  the  BSD-2000 


^A  similar  movable  multiapplicator  phased  array  exists  [12]  and  presumably  could  be  modified  to 
have  adaptive  capability. 
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Figure  1.  A  clinical  hyperthermia  phaaed-array  antenna  system  (photograph  courtesy  of 
BSD  Medical  Corporation). 
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system  are  short  dipoles  with  semiconductor  diode  detectors  [35-37].  A  potentially  better  approach 
to  performing  adaptive  hyperthermia  is  to  use  a  channel-covariance-matrix-based  algorithm  and 
suitable  transmit  calibration  to  control  the  weights.  If  the  electric-field  phase  and  amplitude  are 
measured  (by  modification  of  the  E-field  probe  hardware  to  form  in-phase  and  quadrature  signals), 
the  sample  matrix  inversion  (SMI)  algorithm  [2]  could  be  used  to  form  the  nulls.  The  concept  of 
adaptive  nulling  as  it  applies  to  hyperthermia  will  now  be  described. 

Recently,  a  technique  called  focused  near-field  adaptive  nulling  [38-44]  has  been  investigated 
for  testing  the  performance  of  radar  and  communications  systems.  A  calibration  pr">be  antenna  is 
used  to  focus  the  main  beam  at  approximately  one  diameter  D  of  the  phased  array.  Additional 
(auxiliary)  probe  eintennas  are  located  in  the  sidelobe  region  of  the  phased-array  quiescent-radiation 
pattern.  Adaptive  nulls  are  formed  in  the  direction  of  these  auxiliary  probes.  Computer  simulations 
show  that  a  radar  system  can  be  tested  by  using  the  focused  near-field  nulling  technique  [43]; 
experimental  measurements  of  focused  near-field  adaptive  nulling  have  been  performed  [44].  For 
the  radar  system,  interference  signals  entering  through  uncontrolled  sidelobes  can  reduce  the  signal- 
to-noise  ratio  (SNR);  adaptive  nulling  [1]  is  commonly  used  to  counteract  this  possible  degradation. 
In  the  hyperthermia  system,  high  transmit-antenna  electric-field  strengths  can  create  hot  spots  in 
the  target  tissue.  These  high-temperature  regions  can  possibly  be  alleviated  with  adaptive  nulling 
of  the  electric  field  at  specified  positions  on  the  surface  of  the  target. 

In  the  proposed  adaptive  hyperthermia  concept,  electric-field  nulls  are  used  to  reduce  the 
power  delivered  to  potential  hot  spots.  Computer  simulations  show  that  noninvasive  electric-field 
probes  can  be  used  to  eliminate  hot  spots  that  are  interior  to  the  phantom  target  tissue  [3-6]. 
For  example,  an  eight-element  ring  array  with  60-cm  diameter  surrounding  a  transverse  cross 
section  at  the  prostate  level  is  depicted  in  Figure  2(a).  For  computer  simulation,  the  heterogeneous 
target  is  replaced  with  a  homogeneous  elliptical  phantom.  The  transmit  amplitude  and  phase  of 
each  element  are  assumed  to  be  adaptively  controlled  so  that  a  deep-seated  tumor  site  receives 
maximum  electric-field  power  while  minimizing  the  electric-field  power  at  four  positions  on  the 
target  surface  as  shown  in  Figure  2(b).  Figure  3  summarizes  the  thermal  simulation  for  a  120-MHz 
RF  hyperthermia  ring  array  with  eight  adaptively  controlled  dipole  elements  transmitting  through 
a  constant-temperature  water  bolus  into  a  homogeneous  elliptical  target  region.  The  tumor  site  is 
assumed  to  be  at  the  center  of  the  ellipse.  The  thermal  distribution  in  Figure  3(a)  contains  two 
hot  spots  to  the  left  and  right  of  the  focus.  Note:  a  similar  thermal  distribution  has  been  reported 
in  the  literature  [12].  In  Figure  3(b),  adaptive  nulling  at  four  independent  positions  is  in  effect  and 
the  hot  spots  are  eliminated. 

Consider  a  potential  clinical  application  of  the  proposed  adaptive  hyperthermia  system  with 
the  BSD-2000  Sigma  60  system.  The  BSD-2000  hyperthermia  system  currently  uses  eight  EP-400 
noninvasive  electric-field  probes  to  monitor  clinical  hyperthermia  treatments.  The  eight  EP-400 
probes  (or  eight  EP-lOO  probes)  could  provide  feedback  signals  to  the  adaptive  algorithm.  In 
theory,  with  four  transmit  channels  three  independent  adaptive  nulls  can  be  formed  using  any 
three  of  the  eight  feedback  signals.  There  are  various  methods  to  select  the  desired  null  sites.  The 
adaptive  algorithm  could  minimize  the  three  largest  measured  electric  field  signals  on  the  target 
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Figure  2.  (a)  Transverse  cross  section  at  the  prostate  level.  An  eighUelement  hgperther~ 

mia  ring  array  of  dipole  elements  is  used  to  irradiate  a  target  body  with  sufficient  power 
to  elevate  the  tumor  temperature  in  the  target  body  to  a  therapeutic  level,  (b)  Elliptical 
phantom  target  used  to  model  the  prostate  (target)-level  cross  section.  Shown  are  four 
noninvasive  auxiliary  E~field  probes  that  are  used  in  forming  null  zones  with  widths  con¬ 
trolled  by  the  individual  null  strengths.  The  null  zones  extend  into  the  body  to  reduce  the 
internal  electric  field  at  specified  locations. 


5 


tTMMO 


-15.0  -7.5  0.0  7.5  15.0 


x(cm) 

(b) 

Figure  3.  Simulated  two-dimensional  thermal  pattern  tn  elliptical  phantom  muscle-tissue 
target  surrounded  with  1(FC  constant-temperature  water  bolus.  The  incident  RF  power 
distribution  is  at  J20  MHz  and  the  initial  temperature  of  the  phantom  is  C.  Temper¬ 
ature  contour  levels  are  given  in  2“  C  steps,  (a)  Before  nulling;  hot  spots  on  the  left  and 
right  sides  of  the  target  are  present,  (b)  After  nulling;  hot  spots  are  eliminated. 
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surface.  If  two  or  three  positions  are  identified  as  potential  hot  spots,  adaptive  nulls  could  be 
automatically  formed  at  those  positions.  Alternately,  if  during  the  heating  therapy  a  patient  can 
localize  a  painful  tissue  area,  the  nearest  electric-field  probe  can  be  used  as  the  feedback  signal  to 
the  adaptive  nulling  algorithm.  In  some  situations  more  than  three  independent  adaptive  nulls  may 
be  needed  to  produce  a  therapeutic  thermal  distribution.  An  extension  of  the  present  four-channel 
BSD-2000  system  to  eight  channels  would  provide  up  to  seven  independent  adaptive  nulls  plus  a 
desired  focus. 

The  remainder  of  this  report  is  organized  as  follows.  Section  2  discusses  the  concept  for 
a  noninvasive  adaptive  hyperthermia  system  and  describes  in  detail  the  gradient-search  feedback 
algorithm  used  in  implementing  the  adaptive  array.  Section  3  describi.  *he  materials  and  methods 
used  in  obtaining  the  measured  results.  The  phantom  targets  measured  in  this  report  are  a  saline- 
filled  cylindrical  phantom  and  a  beef  phantom.  Section  4  gives  measured  results  for  an  adaptive 
four-channel  hyperthermia  ring  array  operating  at  CW  frequencies  of  100  MHz  and  120  MHz.  The 
measured  received  RF-power  distributions  before  and  after  adaptive  nulling  with  short-dipole  field 
probes  are  presented.  Section  5  presents  the  conclusions. 
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2.  THEORY 


2.1  Noninvasive  Adaptive  Hyperthermia  System  Concept 

The  concept  of  a  noninvasive  adaptive-nulling  hyperthermia  system  [3-6]  is  shown  in  Figure  4. 
Theoretically,  to  generate  the  desired  field  distribution  in  a  clinical  adaptive  hyperthermia  system, 
receiving  sensors  are  positioned  as  close  as  possible  to  the  focus  (tumor  site)  and  where  high 
temperatures  are  to  be  avoided  (such  as  near  the  spinal  cord  and  scar  tissue).  A  noninvasive 
adaptive  nulling  system  is  achieved  by  placing  auxiliary  sensors  1, 2,  ■  ■  ■ ,  Naux  on  the  target  skin 
as  shown.  The  null  zones  centered  at  each  auxiliary  probe  naturally  extend  into  the  elliptical 
target  region  to  eliminate  undesired  hot  spots.  The  width  of  each  null  zone  is  directly  related 
to  the  strength  of  each  null.  The  strength  of  each  null  (sometimes  referred  to  as  the  amount  of 
cancellation)  is  directly  related  to  the  SNR  at  the  sensor  position.  A  low  SNR  produces  a  small 
amount  of  nulling,  a  high  SNR  a  large  amount  of  nulling.  The  resolution  or  minimum  spacing 
between  the  focus  and  null  positions  is  normally  equal  to  the  half-power  beamwidth  of  the  antenna. 
The  resolution  is  enhanced  somewhat  by  using  weak  nulls  whenever  the  separation  between  the 
null  and  focus  is  closer  than  the  half-power  beamwidth.  The  half-power  angular  beamwidth  of  an 
antenna  aperture  with  diameter  D  in  wavelengths  is  approximated  by 

^HPBW  =  ^ 

where  A  is  the  wavelength.  The  antenna  half-power  focal  beamwidth  (spot  size)  in  units  of  length 
is  expressed  as 


« =  ^HPBW  X  R  ,  (2) 

where  R  is  the  focal  distance  of  the  antenna.  Using  Equation  (1)  and  substituting  R  =  D/2  for  a 
ring  array  focused  at  the  origin  in  Equation  (2)  yields 


Thus,  the  focal  spot  size  or  resolution  of  a  ring  array  is  one-half  the  wavelength  in  the  target  body. 

Initially,  the  hyperthermia  array  is  focused  to  produce  the  required  held  intensity  at  the 
tumor.  An  invasive  probe  is  required  to  achieve  the  optimum  focus  at  depth.  To  avoid  hot  spots,  it 
is  necessary  to  minimize  the  power  received  at  the  desired  null  positions  and  to  constrain  the  array 
weights  to  deliver  a  required  amount  of  transmitted  or  focal-region  power.  The  adaptive  array 
weights  (with  gain  g  and  phase  ^j>)  are  controlled  by  either  the  SMI  algorithm  or  a  gradient-search 
algorithm  to  rapidly  form  the  nulls  before  a  significant  amount  of  target  heating  takes  place.  With 
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Figvn  4-  Noninvasive  adaptive  kyperikermia  system  concept. 


this  adaptive  technique,  it  should  be  possible  to  avoid  hot  spots  and  maintain  a  therapeutic  thermal 
dose  distribution  at  the  tumor. 

2.2  Adaptive  IVansmit- Array  Formulation 

Consider  a  hyperthermia  array  with  N  identical  antenna  elements.  Note:  N  will  represent 
N  pairs  of  identical  antenna  elements  for  the  BSD-2000  system.  The  input  signal  to  each  of  the 
N  array  elements  is  obtained  from  the  weighted  signal  distributed  by  a  power  divider  network. 
The  number  of  adaptive  channels  is  denoted  N.  Let  w  =  (wi,W2,  ■  •  ■  denote  the  adaptive 

channel  weight  vector  as  shown  in  Figure  4.  (Superscript  T  means  transpose.) 

For  an  adaptive  annular  array  focused  at  the  origin  in  homogeneous  tissue,  the  normalized 
quiescent  weight  vector  is  simply  to,  =  (1, 1, 1,  •  •  • ,  1)^;  that  is,  the  amplitude  and  phase  illumi¬ 
nation  are  uniform.  Commonly,  the  weight  vector  is  constrained  to  deliver  a  required  amount  of 
power  to  the  hyperthermia  array  or  to  the  tumor.  For  simplicity  in  the  experimental  adaptive- 
hyperthermia-array  control  software,  the  weights  are  constrained  such  that 

EKNa:  ,  (4) 

where  is  the  transmit- weight  amplitude  for  the  nth  adaptive  channel  and  A"  is  a  constant.  To 

generate  adaptive  nulls,  the  transmit  weights  (phase  and  gain)  are  controlled  by  either  the  SMI 
algorithm  or  a  gradient-search  algorithm.  The  SMI  algorithm  has  the  flexibility  to  operate  in  either 
open-  or  closed-loop  feedback  modes  [45];  the  gradient-search  algorithm  operates  only  in  a  feedback 
mode. 


2.2.1  Gradient- Search  Algorithm 

Gradient-search  algorithms  are  commonly  used  in  adaptive-array  applications  where  the  chan¬ 
nel  correlation  cannot  be  calculated  or  measured.  With  a  gradient  search,  only  the  output  power 
of  the  receiver  channels  needs  to  be  measured  and  is  used  as  a  feedback  signal  to  the  algorithm.  A 
wide  variety  of  gradient  searches  exist  [46-50]. 

Under  conditions  where  only  the  probe-received  power  is  meetsured  (as  in  the  BSD-2000 
system),  it  is  appropriate  to  consider  a  gradient-search  algorithm  to  minimize  the  Bi-field  power  at 
selected  positions.  The  gradient  search  is  used  to  control  the  transmit  weights  iteratively  so  that  the 
RF  signal  received  by  the  probe  array  is  minimized.  The  transmit-array  weights  (gain  and  phase) 
are  adaptively  changed  in  small  increments  2md  the  probe-array  output  power  is  monitored  to 
determine  weight  settings  that  reduce  the  output  power  most  rapidly  to  a  null.  The  mathematical 
formulation  for  the  gradient  search  is  developed  in  a  straightforward  manner  [50,  51]  and  will  now 
be  described  in  the  context  of  hyperthermia.  Although  the  mathematical  formulation  is  given  as  a 
minimization  problem,  the  equations  are  readily  converted  to  the  maxinndzation  problem. 
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The  summation  of  power  received  at  the  electric-field  probes  is  denoted  by  The  adaptive 
array  cancellation  ratio,  denoted  C,  is  defined  here  as  the  ratio  of  the  summation  of  probe-received 
power  after  adaption  Pa  to  the  sununation  of  probe-received  power  before  adaption  pj;  that  is, 

C=—  .  (5) 

Pb 

Consider  now  J  sets  (or  iterations)  of  N  transmit  weights  that  are  applied  to  an  adaptive 
hyperthermia  phased-array  antenna.  In  terms  of  adaptive  nulling,  the  optimum  transmit-weight 
settings  (from  the  collection  of  J  sets  of  N  transmit  weights)  occur  when  the  total  interference  power 
received  by  the  auxiliary  probe  array,  denoted  p’'®'^,  is  minimized.  For  notational  convenience  let  a 
figure  of  merit  F  denote  and  employ  a  method  of  steepest-descent  gradient  search  to  find  the 
optimum  transmit  weights  to  minimize  F;  that  is, 

Fgpt  =  min(Fj)  j  =  1,2,  J  (6) 


Assume  that  there  are  N  complex  transmit  weights  in  the  hyperthermia  phased  array  as 
suggested  by  the  amplitude  and  phase  scatter  diagram  depicted  in  Figure  5.  The  nth  transmit 
weight  in  the  jth  configuration  (or  iteration)  of  transmit  weights  is  denoted 

tt>nj  =  ,  (7) 

where  A„j  is  the  transmit-weight  amplitude  distributed  over  the  range  Amin  to  Amax  and  ^nj  is 
the  transmit-weight  phase  distributed  over  the  range  to  ^max-  The  goal  is  to  find  the  values 
of  amplitude  and  phase  for  each  of  the  N  transmit  weights  such  that  the  figure  of  merit  (p’’®‘^)  is 
minimized.  When  the  figure  of  merit  is  minimized,  adaptive  radiation  pattern  nulls  will  be  formed 
at  the  auxiliary  sensor  positions. 

Assuming  an  initial  setting  of  the  N  transmit  weights,  the  weights  are  adjusted  by  dithering 
them  until  the  optimum  figure  of  merit  is  achieved.  The  goal  is  to  find  the  collective  search 
directions  for  the  N  transmit  weights  such  that  F  decreases  most  rapidly;  that  is,  transmit  weights 
are  selected  so  that  the  directional  derivative  is  minimized  at  (Aj,^j),  where  Aj  and  are  the 
amplitude  and  phase  column  vectors,  respectively. 

The  directional  derivative  of  Fj  is  expressed  in  terms  of  the  amplitude  and  phase  changes  of 
the  transmit  weights  as 


^  ap.  ap. 

n=l  O'Vnj 


(8) 
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Figure  5.  Scatter  diagram  in  rectangular  coordinates  for  the  amplitude  and  phase  of  the 
transmit  weights  in  a  hgperthermia  phased-arrag  antenna.  The  nth  transmit  weight  in  the 
jth  configuration  of  transmit  weights  is  denoted  w„j. 


where  d  means  partial  derivative  and  rjinj,r^nj  are  the  (^,^)  directions  for  which  Fj  is  decreasing 
most  rapidly.  The  directions  r^nj ,  r*„j  are  constrained  by 


N 


+  4nj)  =  1 


(9) 


The  goal  is  to  minimize  D{Fj)  subject  to  the  above  constraint  equation. 

Using  Lagrange  multipliers,  construct  the  Lagrangian  function 

n=l  n=l 

where  G  is  a  constant  to  be  determined.  The  requirement  that  Lj  be  an  extremum  implies 


dLj 

drAnj 

II 

^GtAnj  —  0  »ri  —  li2,  • 

■■,N 

(11) 

dLj 

II 

-  2Gr*„j  =  0  ,n  =  1,2,- 

■■,N  , 

(12) 

or 


1  dFj 

^Anj  -  20dAnj 


1  dFj 
“  2Gd*ni 


Squaring  Equations  (13)  and  (14)  and  invoking  Equation  (9)  yields 


N 


N 


j:aj + -}n/) = 1 = 4^  + (#:)’) 


nsl 


nsl 


^dA 


nj 


nj 


(13) 

(14) 


(15) 


N 


nsl 


^dAn/  ' 


Substituting  Equation  (16)  into  Equations  (13)  and  (14)  gives 


(16) 


(17) 


(18) 


In  Equations  (17)  and  (18)  the  minus  sign  was  chosen  corresponding  to  the  direction  of  maximum 
function  decrease.  (Note:  by  changing  the  minus  sign  to  a  plus  sign  in  Equations  (17)  and  (18), 
the  search  directions  then  correspond  to  the  direction  of  maximum  function  increase,  i.e.,  the  plus 
sign  is  used  to  maximize  the  power  delivered  to  the  focus  or  tumor  site.)  The  partial  derivatives 


dFj  dFj 
dAnj 


;n=l,2,--,lV 


(19) 


represent  the  gradient  directions  for  maximum  function  decrease. 

Because  the  figure  of  merit  F  is  measured  and  cannot  be  expressed  in  an2tlytical  form,  the 
partial  derivatives  are  numerically  evaluated  using  finite  differences.  Thus, 


dFj  ^  AF^nj 

dAnj  2Ai4ny 

dFj  _  AF^nj 
d^nj 

where,  as  shown  in  Figure  6, 

AF^nj  —  ^ji-^nj  "t"  AAnj'j^nj) 


(20) 

(21) 

ji-^nj  ~  AAnj]^nj)  (22) 


AF^j  —  Fj(^Anj',^nj  4"  ^^nj)  ~  f j(,^nj]^nj  ^^nj) 


(23) 
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and  AAnj  and  A^nj  are  the  maximum  step  sizes.  Assume  for  now  that  the  increments  AAnj  and 
A*„j  depend  on  the  iteration  number  j  and  transmit  element  index  n.  Substituting  Equations  (20) 
and  (21)  into  Equations  (17)  and  (18)  gives  the  desired  result  for  the  search  directions 


^Anj  —  ~ 


(24) 


r*nj  =  - 


(25) 


The  new  amplitude  and  phase  settings  of  the  (j  +  l)th  transmit-weight  configuration  are  computed 
according  to 


■4nj+l  =  Anj  +  AAnjrAnj  (26) 

=  ^nj  +  A^njr^nj  j27) 

For  the  current  software  implementation  of  the  gradient  search  in  these  experiments,  assume  (for 
convenience)  that  the  step  sizes  are  independent  of  both  the  iteration  number  and  the  adaptive 
channel  number;  that  is, 


AA„j  =  AA 


(28) 


A^nj  =  A$ 


(29) 


In  some  situations  it  may  be  desirable  to  change  the  step  size  at  each  iteration  as  considered  in 
Farina  and  Flam  [46],  but  that  possibility  has  not  been  explored  in  these  measurements. 

2.3  System  Considerations 

Figure  7  shows  a  block  diagram  for  an  adaptive  hyperthermia  system  controlled  by  the 
gradient-search  algorithm.  The  transmit  weights  wij,- ••  ,Wnjr  ■  •  ,WNj  at  the  jth  iteration  are 
shown  at  the  top  of  the  figure.  The  transmit  phased-array  antenna  induces  a  voltage  across  the  ter- 
minak  of  the  ith  receive  field  probe  antenna.  For  any  given  configuration  of  the  transmit  weights, 
each  weight  is  dithered  by  a  small  amount  in  amplitude  and  phase  and  the  received  powers  at  the 
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AMPLITUDE 


Figure  6.  Figure  of  merit  tuiih  transmiUweight  dithering  for  optimum  search  directions. 
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Figurt  7.  GradieuUaearch  algorithm  block  diagram  for  adapiive~nalling  hgperihermia 
sgstem. 


18 


electric-field  probes  are  stored  in  the  computer  for  calculation  of  the  figure  of  merit,  search  direc¬ 
tions,  and  updated  (j  +  l)th  transmit- weight  configuration.  The  weight  dithering  of  one  transmit 
weight  must  be  done  with  the  remaining  transmit  weights  in  their  jth  state.  The  figure  of  merit 
Fj  in  the  adaptive  hyperthermia  system  is  the  power  received  by  the  auxiliary  probe  array,  as  indi¬ 
cated  in  the  block  diagram.  The  figure  of  merit  is  a  rectangular  matrix  with  dimensions  ( A'  x  4).  A 
dimensionality  of  four  is  due  to  the  plus  and  minus  dithering  of  both  amplitude  and  phase.  Search 
directions  for  the  adaptive  transmit  weights  are  based  on  minimizing  the  auxiliary-probe-array 
received  power  and  are  computed  based  on  Equations  (24)  and  (25).  The  transmit  weights  for 
the  next  configuration  (j  -t-  1)  are  computed  from  Equations  (26)  and  (27).  The  adaptive  weight 
vector  ii>a  is  achieved  when  the  (j  -t-  l)th  weight  configuration  has  converged.  Based  on  computer 
simulations  not  shown  here,  convergence  occurs  in  less  than  one  hundred  iterations,  depending  on 
the  step  sizes  (AA  and  A^)  and  number  of  adaptive  null  positions  A'aux- 
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3.  MATERIALS  AND  METHODS 


The  hyperthermia  phased-array  system  used  in  these  measurements  is  the  BSD  Medical  Cor¬ 
poration  Model  BSD-2000  with  Sigma-OO  applicator.  The  applicator  has  a  59-cm  array  diameter 
with  eight  uniformly  spaced  dipole  elements  (dipole  length  is  44  cm)  operating  over  the  frequency 
band  60-120  MHz  [34],  The  eight  dipoles  are  fed  as  four  active  pairs  of  elements  with  up  to 
500  W  average  power  per  channel.  Each  of  the  four  active  channels  has  an  electronically  controUed 
variable-phase  shifter  for  phase  focusing  the  array.  The  variable  transmit  gain  and  phase  modiiles 
are  controlled  by  digital  to  analog  (D/A)  converters.  Temperature  and  electric-field  probe  sensors 
(both  invasive  and  noninvasive)  are  used  to  monitor  the  treatment.  The  bolus  is  filled  with  cir¬ 
culating  deionized  water,  which  has  a  very  low  RF  propagation  loss.  A  photograph  of  the  SUNY 
Health  Science  Center  Hyp>erthermia  IVeatment  Facility  is  shown  in  Figure  8.  The  hyp>erthermia 
array  is  controlled  and  monitored  by  a  PDOS-based  computer  system  and  the  array  applicator  is 
located  in  an  RF-shielded  room.  Figure  9  shows  the  Sigmar60  ring-array  apphcator,  which  can  be 
moved  by  hand  to  the  desired  longitudinal  position  on  the  patient-support  assembly.  A  photograph 
of  the  four-channel  RF  transmitter  rack  is  shown  in  Figure  10.  Figure  11  shows  the  coaxial  RF 
transmit  feed  lines  for  the  dipole  array. 

Three  models  of  electric-field  probes  are  used  in  these  measurements.  For  invasive  electric-field 
measurements,  a  BSD  Model  EP-500  probe  is  used.  For  noninvasive  electric-field  measurements, 
BSD  Model  Numbers  EP-lOO  and  EP-400  are  used.  The  design  of  the  EP-500  probe  is  as  follows 
[37]:  The  short  metallic  leads  of  a  Schottky  diode  are  connected  between  two  highly  resistive 
conductive  leads.  A  1-cm  length  of  conducting  wire  is  added  to  one  of  the  resistive  leads  and 
responds  to  linear  i>olarization.  The  leads  of  the  diode  and  the  additional  conducting  wire  act  as 
an  o&et-feed  dipole  antenna.  The  lengths  of  the  dipole  antennas  used  in  the  above  three  probes 
are  as  follows:  EP-500  is  a  1-cm  dipole,  EP-lOO  is  a  3-cm  dipole,  and  EP-400  is  a  3-element  array 
of  dipoles  with  a  total  length  of  23  cm.  Two  Bowman  temperature  probes  (provided  with  the  BSD- 
2000  system)  are  used  to  monitor  the  beef  phantom  temp>eratures.  In  this  report,  assume  that  an 
EP-lOO  probe  will  provide  localized  adaptive  nulls;  thus,  the  EP-lOO  probe  is  used  for  providing 
the  feedback  signals.  ^ 

The  saline-filled  (0.9%  NaCl  in  deionized  water)  cylindrical  phantom  (28-cm  diameter  and 
40-cm  length),  located  on  the  patient-supp>ort  sling,  is  shown  in  Figure  12.  Note:  the  cylindrical 
phantom  is  a  6-gallon  polyethylene  bottle.  In  this  figure,  a  dipole  probe  is  mounted  on  the  surface 
of  the  phantom  and  is  connected  to  the  temperatmre/electric-field  monitoring  unit. 

The  beef  phantom  used  in  these  experiments,  depicted  in  Figiu-e  13,  is  a  tapered  cut  obtained 
from  the  hind  leg  just  above  the  knee.  The  front  face  of  the  beef  has  a  horizontal  width  of  38  cm 


^The  EP-400  field  probes  have  also  been  used  for  feedback  signals  and  the  measured  null  strengths 
are  similar  to  the  EP-lOO  measured  data  presented  in  this  report. 
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Figure  8.  Hyperthermia  treatment  facility  at  the  Health  Science  Center,  State  University 
of  New  York. 
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Figure  9.  BSD-2000  Sigma-60  applicator  on  patient-support  assembly  at  SUNY. 
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Figure  11.  BSD-2000  Sigma-60  applicator  with  coaxial  feedlines. 


25 


Figure  12.  Saline-filled  cylindrical  phantom  with  dipole  probe  and  temperature/electric- 
field  monitoring  unit. 


and  a  vertical  height  of  23  cm.  The  thickness  of  the  beef  is  15  cm  and  the  weight  is  10.9  kg.  The 
probe  positions  used  in  the  experiments  are  indicated. 

The  estimated  dielectric  properties  of  the  various  materials  used  in  these  experiments  are 
given  in  Table  1. 


TABLE  1 

Dielectric  Properties  of  Materials  Used  in  the  Experiments 


Material 

Relative  Dielectric  Constant 

Electrical  Conductivity 

Deionized  Water 

78.0 

0.0001 

0.9  %  NaCI  Phantom 

75.0 

1.5 

Polyethylene 

2.25 

0.0002 

Beef  Tissue 

50.0 

l.C 

Fat 

7.5 

0.067 

Bone 

75 

0.067 

A  computer  program  that  implements  the  gradient-search  algorithm  was  developed  by  MIT 
Lincoln  Laboratory  and  integrated  with  the  BSD  Medical  Corporation  quad-amplifier  control  and 
temperature/electric-field  monitoring  software  at  SUNY.  The  MIT/SUNY  integrated  software  al¬ 
lows  the  user  to  run  either  adaptive-nulling  or  adaptive-focusing  algorithms  and  to  monitor  the 
output  power  of  the  electric-field  probes  during  the  gradient  search.  The  softwsu'e  also  writes  mea¬ 
sured  data  to  a  file.  The  adaptive-nulling  algorithm  uses  phase  and  amplitude  control  to  form  the 
nulls — the  adaptive-focusing  algorithm  uses  only  phase  control  to  mtiximize  the  focal-point  electric 
field. 
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•  ELECTfUC-FIELD  PROBE  SneS 
X  TEUPERAWRE  PROBE  SrrES 

FRONT  PACE  VIEW 


THICta^SSm  15cm 


Figure  IS.  Front  face  dimensions  and  probe  positions  for  beef  phantom. 


4.  MEASURED  RESULTS 


In  this  section,  measured  data  are  presented  for  the  saline-iilled  cylindrical-phantom  and 
beef-phantom  targets.  Measurement  frequencies  were  100  MHz  for  the  cylindrical  phantom  and 
120  MHz  for  the  beef  phantom.  These  frequencies  were  selected  by  monitoring  the  reflected  power 
level  of  the  transmitter  channels  as  a  function  of  frequency.  Frequencies  where  the  reflected  power 
is  low  were  determined  by  trial  and  error. 

The  experimental  test  setup  for  the  saline-filled  cylindrical  phantom  is  shown  in  Figure  14. 
Three  electric-field  probes  are  used  to  monitor  the  electric-field  distribution.  One  probe  is  mounted 
invasively  at  the  center  of  the  cylindrical  phantom  and  two  probes  are  mounted  on  the  left  and 
right  sides  of  the  phantom.  The  power  measured  at  the  invasive  probe  site  is  taken  as  the  power 
delivered  to  a  fictitious  tumor  site. 

The  experimental  test  setup  for  the  beef  phantom  is  shown  in  Figure  15.  Three  electric-field 
probes  are  used  to  monitor  the  electric-field  distribution.  One  probe  is  mounted  invasively  at  the 
center  of  the  beef  phantom  and  two  probes  are  mounted  on  the  left  and  right  sides  of  the  phantom. 
The  power  measured  at  the  invasive  probe  site  is  taken  as  the  power  delivered  to  a  fictitious  tumor 
site. 

For  the  first  experiment,  one  adaptive  null  is  formed  on  the  right  side  of  the  cylindrical 
phantom  while  the  electric  field  is  monitored  at  the  center  and  left  side  of  the  phantom.  In  the 
second  experiment,  two  adaptive  nulls  are  formed  independently  on  the  left  and  right  sides  of  the 
cylindrical  phantom  while  monitoring  the  tumor-site  power.  The  longitudinal  characteristics  of  one 
of  the  adaptive  nulls  is  quantified  by  measurement.  In  the  third  experiment,  an  adaptive  focusing 
gr2ulient  search  was  run  for  the  cylindrical  phantom.  The  focus  was  positioned  on  the  surface  of 
the  phantom.  In  the  fourth  experiment,  one  adaptive  null  on  the  side  of  the  beef  is  measured. 
Thermocouple  data  are  presented  for  this  case. 

The  appendix  lists  the  raw  A/D  data  for  the  measured  electric-field  probe  powers  in  the 
phantom  experiments.  In  the  following  sections,  the  A/D  converter  data  is  converted  to  power  in 
dB  by  computing  lOlog(p)  where  p  is  the  value  (meeisured  power)  of  the  A/D  converter. 

4.1  One  Adaptive  Null  for  Cylindrical  Phantom 

The  invasive  EP-500  probe  is  mounted  in  the  center  of  the  cylindrical  phantom  and  measures 
the  electric-field  power  at  the  simulated  tumor  site.  The  probe  is  inserted  in  a  catheter  through 
the  rubber  plug  on  top  of  the  phantom.  An  EP-400  probe  is  taped  to  the  left  side  of  the  phantom 
and  an  EP-lOO  probe  is  taped  to  the  right  side  of  the  phantom.  The  EP-400  probe  measures 
a  reference  electric  field,  while  the  EP-lOO  measures  the  adaptive-nulling  feedback  signal.  A  50- 
iteration  adaptive-nulling  gradient  search  was  executed  for  this  test  configuration  at  100  MHz.  The 
gradient-search  step  sizes  in  terms  of  raw  D/A  converter  states  are  AA  =10  and  =50.  Note:  the 
D/A  converters  have  12  bits,  thus  4096  states.  The  commanded  initial  transmit-weight  amplitude 
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Figure  14-  Experimental  test  setup  for  the  saline-filled  cylindrical  phantom.  Three  probes 
are  used  to  monitor  the  electric  field  during  the  adaptive-nulling  measurements. 
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Figure  15.  Beef  phantom  target  with  electric-field  probes  and  thermocouples. 
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states  on  channels  1  to  4  ^  are  (commanded  Li/A  states:  500,  500,  500,  500)  or  (forward  power  in 
watts:  215  W,  215  W,  215  W,  215  W)  “*  and  the  corresponding  phase  states  are  (commanded  D/A 
states:  500,  500,  500,  500)  or  (equivalent  phase  lag;  89°,  89°,  89°,  89°).  After  50  iterations  of  the 
nulling  algorithm,  the  transmit-amplitude  states  are  (commanded  D/A  states:  620,  668,  417,  294) 
or  (forward  power  in  watts:  370  W,  440  W,  98  W,  8  W)  and  the  phase  states  are  (commanded  D/A 
states:  58,  578,  594,  664)  or  (equivalent  phase  lag:  45°,  94°,  95°,  99°).  The  measured  electric  field 
as  a  function  of  iteration  number  is  displayed  in  Figure  16.  There  are  three  curves  plotted;  each 
curve  is  normalized  to  0  dB  at  the  initial  iteration.  The  null-site  probe  power  clearly  decreases 
during  the  first  40  iterations.  After  44  iterations,  the  average  null-site  amplitude  (cancellation) 
is  about  —20  dB.  The  tumor  site  power  has  dropped  by  about  5  dB  at  iteration  50;  this  power 
reduction  is  attributed  in  part  to  the  beam  peak  shifting  away  from  the  focus.  The  reference-site 
power  stays  relatively  constant  during  the  50  iterations. 

4.2  Two  Adaptive  Nulls  for  Cylindrical  Phantom 

The  second  experiment  involves  forming  two  simultaneous  adaptive  nulls.  The  invasive  EP- 
500  probe  is  mounted  in  the  center  of  the  cylindrical  phantom  and  measures  the  electric-field  power 
at  the  simulated  tumor  site.  Two  EP-lOO  probes  are  taped  to  the  left  and  right  sides  of  the  phantom 
and  measure  the  adaptive-nulling  feedback  signals.  A  50-iteration  adaptive-nulling  gradient  search 
was  executed  for  this  test  configuration  at  100  MHz.  As  before,  the  gradient-search  step  sizes  in 
terms  of  raw  D/A  converter  states  are  AA  =10  and  A$  =50.  The  initial  transmit-weight  amplitude 
states  on  channels  1  to  4  are  (commanded  D/A  states:  500,  500,  500,  500)  or  (forward  power  in 
watts;  215  W,  215  W,  215  W,  215  W)  and  the  corresponding  phase  states  are  (commanded  D/A 
states;  500,  500,  500,  500)  or  (equivalent  phase  lag:  89°,  89°,  89°,  89°).  After  50  iterations  of 
the  nulling  adgorithm,  the  transmit-amplitude  states  are  (commanded  D/A  states:  744,  379,  482, 
396)  or  (forward  power  in  watts;  510  W,  55  W,  190  W,  76  W)  and  the  corresponding  phase  states 
are  (commanded  D/A  states:  1,  676,  228,  781)  or  (equivalent  ph^lse  lag:  0°,  99°,  58°,  101°).  The 
measured  electric  field  as  a  function  of  iteration  number  is  displayed  in  Figure  17.  The  two  null-site 
probe  powers  decrease  nearly  monitonically  over  the  50  iterations.  At  iteration  50,  the  cancellation 
for  null-site  probe  #2  is  -18.0  dB,  while  the  cancellation  for  null-site  probe  #3  is  —11.5  dB.  The 
tumor-site  probe  power  stays  relatively  constant  during  the  50  iterations.  With  the  two  symmetric 
nulls,  the  main  beam  tends  to  stay  centered  at  the  desired  focus  [3]. 


®The  numbering  of  transmit-array  elements  1, 2,3,4  correspond  to  the  anterior,  posterior,  left,  right 
quadrants,  respectively,  depicted  in  Figure  16.  Note  that  "right  side”  refers  the  side  of  the  Sigma-60 
that  contains  the  coaxial  feed  cables. 

‘*The  D/A  states  are  converted  to  power  by  a  lookup  table  constructed  from  calibrated  power-meter 
measurements  performed  at  SUNY. 


32 


RUN  36. 3  DEC  61.  BSO-2000  SYSTEM  AT  SIMY 


ITERATION  NUMBER 


EXPERMENTAL  CONnOURATKM 

ANTERtOR  SIQMA-60 


Figvn  16.  Experimental  demonstration  of  an  adaptive-nulling  hyperthermia  system  at 
SUNY.  The  figure  shows  the  measured  electric  field  as  a  function  of  adaptive-nulling 
gradient-search  iteration  number  for  a  saline-filled  cylindrical  phantom.  One  adaptive 
null  is  formed  at  probe  position  2. 
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Figure  17.  Measured  electric  field  as  a  function  of  adaptive-nulling  gradient-search  iter¬ 
ation  number  for  a  saline-filled  cglindrical  phantom.  Two  adaptive  nulls  are  formed  at 
probe  positions  2  and  S.  Notice  that  the  power  measured  by  the  invasive  electric-field  probe 
is  relatively  constant  over  the  50  iterations. 
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Using  the  adaptive  weights  from  iteration  50,  a  measurement  of  the  electric  field  in  the  vicinity 
of  one  of  the  nulls  was  performed.  Probe  3  was  moved  uniformly  to  nine  positions  longitudinally 
on  the  surface  of  the  cylindrical  phantom;  the  nine  positions  scanned  20  cm.  The  electric  field 
before  and  after  nulling  is  presented  in  Figure  18.  Before  nulling,  the  radiation  pattern  has  broad 
coverage  and  the  half-power  beamwidth  is  at  least  20  cm  or  more.  After  nulling,  the  electric  field  is 
reduced  at  each  measurement  position  over  the  20-cm  distance.  The  electric  field  tends  to  increase 
as  the  probe  moves  away  from  the  null  position  at  y  =  0.  This  measured  result  is  consistent  with 
computer  simulations  where  a  finite- width  null  is  observed  [52]. 

4.3  Adaptive  Focus  for  Cylindrical  Phantom 

The  third  experiment  involves  forming  an  adaptive  focus  on  the  surface  of  the  cylindrical 
phantom.  One  EP-lOO  probe  was  taped  to  the  right  side  of  the  phantom  and  measured  the  adaptive- 
focusing  feedback  signal.  ®  A  30-iteration  adaptive-focusing  gradient  search  was  executed  for  this 
test  configuration  at  100  MHz.  The  gradient-search  step  sizes  in  terms  of  raw  D/A  converter  states 
are  AA  =0  and  =50.  Because  the  amplitude  step  size  is  zero,  only  phase  focusing  is  in  effect. 
The  initial  transmit-weight  amplitude  states  on  channels  1  to  4  are  (comm'aded  D/A  states:  500, 
500,  500,  500)  or  (forward  power  in  watts:  215  W,  215  W,  215  W,  215  W)  and  the  corresponding 
phase  states  are  (commanded  D/A  states:  700,  700,  700,  700)  or  ^equivalent  phase  lag:  100°, 
100°,  100°,  100°).  After  30  iterations  of  the  focusing  algorithm,  the  transmit-amplitude  states  are 
(commanded  D/A  states:  500,  500,  500,  500)  and  the  corresponding  phase  states  are  (commanded 
D/A  states:  1024,  659,  374,  132)  or  (equivalent  phase  lag:  113°,  98°,  78°,  30°).  The  measured 
electric  field  as  a  function  of  iteration  number  is  displayed  in  Figure  19.  The  focal-probe  power 
increases  by  0.8  dB  over  the  30  iterations.  The  data  indicate  convergence  in  about  15  iterations. 

4.4  One  Adaptive  Null  for  Beef  Phantom 

The  invasive  EP-500  probe  is  mounted  in  the  center  of  the  beef  phantom  and  measures  the 
electric-field  power  at  the  simulated  tumor  site.  An  EP-400  probe  is  taped  to  the  left  side  of  the 
phantom  and  an  EP-lOO  probe  is  taped  to  the  right  side  of  the  phantom.  Note:  The  probe  positions 
were  already  shown  in  Figure  13.  The  EP-400  probe  measures  a  reference  electric  field,  whih'  the 
EP-lOO  measures  the  adaptive-nulling  feedback  signal.  A  50-iteration  adaptive-nulling  gradient 
search  was  executed  for  this  test  configuration  at  120  MHz.  The  gradient-search  step  sizes  in  terms 
of  raw  D/A  converter  states  are  AA  =10  and  A$  =50.  The  initial  transmit-weight  amplitude  states 
on  channels  1  to  4  are  (commanded  D/A  states:  500,  500,  500,  500)  and  the  corresponding  pha.se 
slates  are  (commanded  D/A  states:  500,  500,  500,  500)  or  (equivalent  phase  lag:  103°,  103°.  103°. 


®Note:  an  attempt  was  made  to  focus  the  array  with  the  EP-500  invasive  probe  at  the  center  of  the 
phantom.  However,  the  signal  level  is  low  (or  the  array  was  perfectly  focused)  so  no  improvement 
in  focusing  was  measured. 
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Figure  18.  Measured  electric  field  tn  a  longitudinal  cut  on  the  surface  of  the  cylindrical 
phantom  before  and  after  nulling. 
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Figure  19.  Measured  electric  field  as  a  function  of  adaptive-phase-focusing  gradient- 
search  iteration  number  for  a  saline-filled  cylindrical  phantom.  The  power  level  at  the 
focus  increases  by  0.8  dB  due  to  the  adaptive  gradient  search. 
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103°).  After  50  iterations  of  the  nulling  algorithm,  the  trainsmit-amplitude  states  are  (conunanded 
D/A  states:  465,  636,  472,  427)  or  (forward  power  in  watts:  163  W,  390  W,  170  W,  112  W) 
and  the  phase  states  are  (commanded  D/A  states:  278,  195,  188,  988)  or  (equivalent  phase  lag: 
76°,  60°,  58°,  131°).  The  measured  electric  field  as  a  function  of  iteration  number  is  displayed  in 
Figure  20.  The  null-site  probe  power  decreases  during  the  first  thirty  iterations.  After  30  iterations, 
the  average  null-site  amplitude  (cancellation)  is  about  —20  dB.  The  tumor-site  power  has  dropped 
by  about  2  dB  at  iteration  50;  this  power  reduction  is  attributed  in  part  to  the  beam  peak  shifting 
away  from  the  focus  as  observed  for  the  cylindrical  phantom.  The  reference  site  power  is  observed 
to  increase  by  about  2  dB  during  the  50  iterations. 

To  record  calibrated  thermocouple  data,  the  RF  was  turned  off  after  each  10  iterations  ^  for 
about  five  minutes.  During  the  period  when  the  RF  is  off,  the  temperature-monitoring  mode  of 
the  BSD-2000  system  is  accessed  and  the  data  recorded.  The  RF  is  then  turned  back  on  with 
the  transmit  weights  set  to  the  previous  iteration.  The  measured  temperatures  at  the  tumor  site 
and  null  site  are  shown  in  Figure  21.  The  initial  temperatures  were  approximately  equal;  after 
80  minutes  there  is  a  4°C  higher  temperature  at  the  tumor  site  compared  to  the  null  site.  Thus, 
the  data  suggest  that  adaptive  nulling  may  be  effective  in  improving  the  thermal  distribution  in 
hyperthermia.  More  detsdled  thermal  measurements  in  three  dimensions  have  yet  to  be  attempted 
in  this  study. 


^The  length  of  time  for  10  iterations  is  approximately  15  minutes  with  the  present  soft¬ 
ware/hardware  configuration.  A  faster  iteration  speed  is  required  for  practical  use  of  the  gradient 
search  in  patient  therapy.  By  optimizing  the  software/hardware  configuration  for  adaptive  feedback 
operation,  a  faster  speed  should  be  possible. 
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Figure  SO.  Measured  electric  field  as  a  function  of  adaptive-nulling  gradient-search  iter¬ 
ation  number  for  a  beef  phantom.  One  adaptive  null  is  formed  at  probe  position  S. 
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Figure  21.  Measured  temperatures  in  the  beef  phantom.  Eighty  minutes  after  the  start 
of  the  experiment,  a  higher  temperature  exists  at  the  tumor  site  compared  to  the  null 


5.  CONCLUSION 


An  adaptive  hyperthermia  system  has  been  experimentally  demonstrated  by  modifying  a 
commercial  hyperthermia  system  at  the  State  University  of  New  York  Health  Science  Center  at 
Syracuse.  The  measurements  shown  in  this  report  demonstrate  that  the  electric-field  radiation  pat¬ 
tern  of  a  hyperthermia  phased  array  can  be  controlled  by  an  adaptive  algorithm.  Both  electric-field 
nulling  (power  minimization)  and  electric-field  focusing  (power  maximization)  at  desired  phantom 
target  positions  are  experimentally  demonstrated.  The  body  of  this  work  consists  of  electric  field 
measurements;  however,  limited  beef-phantom  temperature  data  were  also  measured  and  shown. 

The  purpose  of  electric-field  nulling  is  to  reduce  the  occurrence  of  undesired  hot  spots  inside 
or  on  the  surface  of  the  target  body.  In  contrast,  electric-field  focusing  is  intended  to  maximize  the 
RF  power  delivered  to  a  tumor  site.  The  data  presented  in  this  report  suggest  that  both  of  these 
goals  may  be  achieved  with  an  adaptive  hyperthermia  phased-array  system. 

In  this  report,  the  adaptive  algorithm  is  a  gradient-search  feedback  technique  (method  of 
steepest  descent /ascent)  derived  in  Section  2.  The  gradient  search  employs  transmit-weight  gain 
and  phase  dithering  together  with  electric-field  probe-power  measurements,  at  the  desired  target 
positions,  to  calculate  the  required  gradient-search  directions  for  RF  nulling.  For  RF  focusing,  a 
phase-only  transmit-weight  dithering  algorithm  was  used. 

ihe  hyperthermia  system  investigated  here  is  the  BSD  Medical  Corporation  Model  BSD- 
2000  with  Sigma-60  phased-array  antenna  applicator.  The  array  consists  of  four  transmit  channels 
connected  to  four  pairs  of  dipole  radiators.  Short-dipole  electric  field  sensors  are  used  to  pro¬ 
vide  feedback  amplitude  signals  for  the  adaptive  array.  Two  phantom  targets  were  investigated 
in  this  report:  a  homogeneous  saline-filled  cylindrical  target  and  a  heterogeneous  (tissue,  bone, 
and  fat)  beef  target.  One  and  two  independent  adaptive  electric-field  nulls  have  been  attempted 
and  demonstrated  with  the  saline-filled  cylindrical  phantom  at  100  MHz.  The  electric-field  longi¬ 
tudinal  surface  distribution  was  measured  before  and  after  adaptive  nulling  with  the  saline-filled 
cylindrical  phantom.  The  longitudinal  field-probe  mettsurements  indicate  a  finite-width  null  char¬ 
acteristic.  An  adaptive  focus  was  measured  with  an  electric-field  probe  positioned  on  the  surface 
of  the  cylindrical  phantom.  One  adaptive  null  was  attempted  and  demonstrated  at  120  MHz  with 
the  beef  phantom.  During  the  adaptive-nulling  beef-phantom  testing,  the  measured  thermocouple 
data  indicate  a  4®C  higher  tumor-site  temperature  compared  to  the  temperature  at  the  surface 
of  the  beef.  The  measured  adaptive  null  strengths  were  on  the  order  of  —10  dB  to  —30  dB.  For 
practical  implementation  of  the  adaptive-nulling  algorithm  the  required  minimum  null  strength  for 
local  near-surface  hot-spot  reduction  is  assumed  to  be  on  the  order  of  —3  dB  to  —6  dB,  which  is 
readily  achieved. 

The  measured  phantom  data  are  encouraging  and  indicate  that  a  potential  improvement  in 
hyperthermia  thermal  dose  exists  with  adaptive  array  control.  Further  experiments  are  planned. 
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